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Abstract. Event-by-event fluctuations of initial QCD-matter density produced in 
heavy-ion collisions at RHIC enhance the production of thermal photons significantly 
in the region 2 < pr < 4 GeV/c compared to a smooth initial-state averaged profile 
in the ideal hydrodynamic calculation. This enhancement is a an early time effect due 
to the presence of hotspots or over-dense regions in the fluctuating initial state. The 
effect of fluctuations is found to be stronger in peripheral than in central collisions. 



1. Introduction 

Photons have long been considered as one of the most powerful probes to characterize 
the initial state of the system produced in ultra-relativistic heavy- ion collisions [1] . This 
is complementary to the pt spectra of bulk hadrons which reflect the late conditions close 
to freeze-out. The thermal emission of photons shows a strong temperature dependence, 
high px photons are mostly emitted from the hot and dense early stage of the system 
evolution when the hydrodynamic flow is weak, whereas low px photons are emitted from 
the flow boosted relatively cold later part of the system. Recent studies have shown 
that event-by-event QCD-matter density fluctuations in the initial conditions (IC) are 
needed (in addition to a correct reference plane definition) for reproducing the measured 
charged-particle elliptic flow even for most central collisions, which was underestimated 
by all earlier hydrodynamic calculations using smooth IC [2]. 

We show that fluctuations in the initial state enhance the production of thermal 
photons significantly for > 1 GeV/c compared to a smooth IC. This makes the 
thermal photons in that pt range especially suitable to probe the hotspots in the IC [3] . 

2. Event-by-event hydrodynamics and thermal photons 

We use the event-by-event hydrodynamical model developed in [2] to calculate the 
production of thermal photons from fluctuating as well as from smooth IC. A Monte 
Carlo Glauber model with the standard two-parameter Woods-Saxon nuclear density 
profile is used to create the initial state nucleon configurations, and to compute the 
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Figure 1. (From [3 . Color online) [Left] Thermal photons from a smooth (long 
dashed blue line) and fluctuating (solid red line) IC along with PHENIX direct photon 
data [H |5] for 0-20% centrality. [Right] a dependence of the results. 

number of wounded nucleons. The initial density profile is taken to be proportional to 
the number of wounded nucleons (WN), where entropy density s is distributed in the 
(x, y) plane around the wounded nucleons using a 2D Gaussian smearing, 

< X >V) =2^5 g ex P( 2^ )• C 1 ) 

Here K is a fixed overall normalization constant and a is a free parameter which 
determines the size of the fluctuations. The initial formation time of the plasma is 
taken as To = 0.17 fm/c, motivated by the EKRT minijet saturation model We 
choose a default value for the size parameter as a = 0.4 fm [2]. The equation of state 
which shows a sharp cross-over transition from the plasma phase to the hadronic phase 
is from [5]. The temperature at freeze-out is taken as 160 MeV (see [2] and [3] for more 
details). This model has been applied successfully to reproduce both the measured 
centrality dependence as well as the pt dependence of charged particle elliptic flow up 
to pt ~ 2 GeV/c, and the pion p T spectra up to ~ 3 GeV/c [2]. 

We use the complete leading order rates R = EdN/d 3 pd 4 x of [6] for the emission 
from plasma and the parametrization given in [7] for the hadronic matter emission 
which at present can be considered as the state of the art for photon production. The 
switching from the plasma rates to hadronic rates is done at a temperature 170 MeV. 
The total thermal emission from the produced QCD matter is obtained by integrating 
the emission rates over the space-time history of the fireball as 

EdN/d 3 p = J d*xR(E*(x),T(x)), (2) 

with T(x) as the local temperature and E*(x) = p^u^x), where p^ is the four- 
momentum of the photon and is the local four-velocity of the flow field. 

3. Results and discussions 

Figure [1] (left panel) shows our results for thermal photon spectra from a smooth and 
from fluctuating IC for 200A GeV Au+Au collisions at RHIC along with PHENIX 
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Figure 2. (From [3]. Color online) Time evolution of thermal photon yield for pt 
values of (a) 1, (b) 3, and (c) 5 GeV/c. Results are compared with an initial state 
averaged event and four different random events from the fluctuating IC. 

data [8j [9] for the 0-20 % centrality bin. The result from the fluctuating IC scenario 
is obtained by averaging photon spectra from a sufficiently large number of random 
events. The smooth profile is obtained by taking an average of 1000 initial profiles 
which is enough to remove essentially all the fluctuations [3]. We see that for both the 
smooth and the fluctuating IC, the spectra are dominated by radiation from the QGP 
phase in the entire p? range shown in the figure. We further notice that even with a 
much smaller freeze-out temperature (~ 120 MeV) we get significant contributions only 
below pt ~ 1.5 GeV/c as shown by the brown dashed dotted line in Figured] 

The slope of the photon px spectrum from the fluctuating IC is about 10% flatter 
compared to the slope from the smooth IC in the region 2 < < 4 GeV/c. For px < 2 
GeV/c, with the fluctuating IC we obtain 20-40% more photons than with the smooth 
profile, whereas for pt > 2 GeV/c the two results differ almost by a factor of 2. We 
know that the thermal emission of photons is exponential in temperature and linear in 
radiating volume. As a result, the hotspots in the fluctuating IC produce more high 
Pt photons than the smooth profile and the difference between the two IC increases 
towards higher values of transverse momentum. We find that the photon results from 
the fluctuating IC show a better agreement with the PHENIX experimental data for 
2 < Pt < 4 GeV/c leaving enough space for the other sources of direct photons besides 
the thermal contribution. 

The thermal photon production from the fluctuating IC is found to be quite 
sensitive to the value of the fluctuation size parameter a, where the enhancement in 
the production compared to a smooth IC is found to be maximal towards smaller values 
of a and also for larger values of pt (see right panel of Figured]). 

The time evolution of thermal photon yield for different values of pt is shown in 
Figure [2] which clearly explains the dynamics leading to the results presented in Figured] 
We see that most of the high pr photons are emitted from very early stage of the system 
expansion and their production drops by a few order of magnitude within a couple of 
fm/c. We further notice that the cold events (or events having less than average entropy) 
produce more high pt photons than the smooth IC for pt > 3 GeV/c. 
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Figure 3. (Color online) [Left] Centrality dependence of the thermal photon spectra 
from smooth and fluctuating IC. [Right] Results at RHIC and LHC are compared. 



The centrality dependence of the results is shown in Figure [3] (left panel) where 
the difference between the smooth and fluctuating IC is found to increase from central 
to peripheral collisions. The size of the produced system becomes smaller towards 
peripheral collisions and the effect of the hot spots in the fluctuating IC becomes more 
pronounced, leading to the shown result. We have also noticed that the effect of initial 
state fluctuations in the excess production of high photons is weaker at LHC than 
at RHIC [10] (see right panel of Figure [3]) . 

In conclusion, we see that fluctuations in the initial density distribution enhance 
significantly the production of thermal photons compared to a smooth initial state 
averaged profile in ideal hydrodynamic calculation. This is an early time effect solely 
due to the presence of hotspots in the fluctuating IC and this effect is found to be 
stronger towards peripheral collisions. This calls for a more detailed study of the effects 
of fluctuations for lower beam energies and for smaller size systems. In addition to that, 
estimation of photon elliptic flow from the fluctuating IC will be very interesting. 
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